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We apply the Green's function based full-potential screened Korringa-Kohn-Rostoker method in 
conjunction with the local density approximation to study the surface energies of the noble and the 
fee transition and sp metals. The orientation dependence of the transition metal surface energies can 
be well described taking into account only the broken bonds between first neighbors, quite analogous 
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to the behavior we recently found for the noble metals [see cond-mat/0105207|. The (111) and (100) 
T-H ' surfaces of the sp metals show a jellium like behavior but for the more open surfaces we find again 

\ the noble metals behavior but with larger deviation from the broken-bond rule compared to the 

^•ZJ . transition metals. Finally we show that the use of the full potential is crucial to obtain accurate 

surface energy anisotropy ratios for the vicinal surfaces. 

-t— > ' 

^ \ PACS numbers: 71.15.Nc, 71.15.Cr, 71.20.Gj 
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I. INTRODUCTION 

^ ! 

^ ' The surface energy is one of the most fundamental solid state properties since it determines the equilibrium shape 
of a mezoscopic crystal and plays a decisive role in phenomena like roughening, faceting and crystal growlh. Despite 
of their importance, surface energies are difficult to determine experimentally and just few data exism. Most of 
these experiments are performed at high temperatures and contain uncertainties of unknown magnitudeEJ. The most 
y—i, comprehensive e:!ajfirimental data stem from surface tension measurements in the liquid form being extrapolated to 
^-H zero temperatureaB, which cannot provide orientation specific information. Gold crystalliteso and single surfacesu have 
attracted a lot of attention aiming to study the orientation dependence of the surface energy, but these experiments, 
as it was also the case for experiments on In and Pb crystallitesQ, are performed at high temperatures so that the 
results are difficult to interpret. Entropy terms, describing the lower vibrational frequencies of the atoms at the 
surface as compared to the bulk, the formation of kinks and finally the creation of holes and pillboxes at the ipw- index 
surfaces, have to be added to the total free energy. At such high temperatures the surface-melting faceting3, i.e. the 
^ break-down of a vicinal surface in a dry and a melted one, plays a predominant role. Also the measureijaent of core 
Q " level shifts at the surface has been proposed as an indirect measurement of the surface energy anisotropyH. Recently, 
O Bonzel and EdmundtsEI have shown that analyzing the equilibrium shape of crystallites at various temperatures by 
^ scanning tunneling microscopy can yield absolute values of the surface energies versus temperature, but this technique 
• , has not yet been applied. 

/\ • During tiieiast years there have been several atteppts to calculate the surface e nejgv of metals using either ab-initio 
^ , techniqueal3"ll3, tight-binding parameterizational^EJ or semi-empirical methodst3lla. Methfessel and collaborators 

- - ' were the first to study rtke trends in the surface energy, work function and relaxation for the whole series of bcc and 
fee 4d transition metalst2l, using a full-potential version of the linear muffin-tin orbitals (LMTO) method in conjunction 
with the local-spin density approximation to the exchangpcorrelation potentialOEa. In the same spirit Skriver and 
co-workers have used a Green 's function LMTO techniqucEj to calculaip the surface energy and the work function of 
most of the elemental metalsOO'EB. Recently, Vitos and coUaboratorsEa using their full-charge density (PCD) Green's 
function LMTO technkuie in the atomic sphere approximation (ASA)c3 in conjunction with the generalized gradient 
approximation (GjGA)o constructed a large database that contains the low-index surface energies for 60 metals in 
the periodic tablec^. Their results present ajmean deviation of 10 % from the full-potential results by Methfessel 
and collaborators forJhe Ad transition metalaHj. Afterwards, they have used this database in conjunction with the 
pair-potential modelEZI to-^timate the formation energy for monoatomic steps on low-index surfaces for an ensemble 
of the bcc and fee metalttl. 

In reference ^ we have demonstrated that the surface energies of noble metals scale accurately with the number of 
broken bonds between first neighbors. This broken-bond rule is very useful for the estimation of the surface energies 
of vicinal surfaces and of the step energies; the latter ones can be calculated as the energy difference between a vicinal 
and a flat surface. In this contribution we investigate the question whether the broken-bond rule can also be applied 
to the surface energies of the other paramagnetic fee metals: the transition metals Rh, Pd, Ir and Pt, and the sp 
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metals Ca, Sr, Al and Pb. To calculate the surface energies we used the recently developed screened Korringa-Kohn- 
Rostoker (KKR) metbed which has been already used to calculate the magnetic properties of Ad monoatomic rows 
on Ag vicinal surfacea23. In Section II we analyze the details of our calculations and the convergence of our results. 
We also discuss the importance of relativistic effects. In Section III, we present the surface energies of the transition 
and sp metals and discuss the applicability of the broken bond rule for these systems. All results presented in section 
III are obtained accounting for relativistic effects in the scalar-relativistic approximation. Finally we discuss the use 
of the full-potential instead of the atomic sphere approximation 

II. METHOD OF CALCULATION 
A. Computational Details 

To perform the calculations, we used the Vosko, Wiik and Nusair parameterizationlil for the local density approx- 
imation (LDA) to the exchange-correlation potentiaO to solve the Kohn-Sham equations within the screened KKR 
method that was recently developed in our groupE3. Its main advantage is that it can treat both 2D and 3D systems 
in the same footing. Both the atomic sphere approximation (ASA)lj and the capability to treat the full-potential 
(FP) are implemented in this scheme. The ASA calculations take into account the full charge density. It was shown 
by Andersen and collaborators that the charge|-d|ensity obtained in this way for spherically symmetric potentials is 
close to the density obtained using a FP methodEd. The full-potential is implemented by using a Voronoi construction 
of Wigner-Seitz polyhedra that fill the space as described in reference ^ A repulsive muffin-tin potential (4 Ry 
highlJs used as reference system to screen the free-space long-range structure constants into exponentially decaying 
onest3. For the screening we took for all metals interactions up to the second neighbors into account leading to a 
tight-binding (TB) cluster around each atom of 19 neighbors. To calculate the charge density, we integtated along a 
contour on the complex energy plane, which extends from the bottom of the band up to the Fermi levelEj. Due to the 
smooth behavior of the Green's functions for complex energies, only few energy points are needed; in our calculations 
we used 22, energy points. For the Brillouin zone (BZ) integration, special points are used as proposed by Monkhorst 
and PackO. Only few tens of k|| are needed to sample the BZ for the complex energies, except for the energies close 
to the real axis near the Fermi level for which a considerably larger number of k|| points is needed. Here we used from 
~ 300 points for the vicinal surfaces up to ~ 800 points for the (110) surface. In addition we used a cut off of imax=Q 
for the mutlipole expansion of the charge density and the potential and a cut off of ^maa;=3 for the wavefuctions. 
Finally in our calculations the core electrons are allowed to relax during the self-consistency. 

To simulate the surface we used a slab with N metal layers and N^^ac vacuum layers from each side. We have 
converged the number of metal and vacuum layers so that our surface energies are converged within 0.01 eV. The 
number of layers needed to converge the surface energies increases with the roughness of the surface and we had to 
use 12 layers of the fee metal for the (111) surface, 14 for the(lOO), 18 for the (110), 21 for the (311), 30 for the (331) 
and the (210) surfaces. We have also used 3 vacuum layers from each side of the slab. For all the systems studied we 
used the experimental lattice parameters: 3.80 A for EJi, 3.89 A for Pd, 3.84 A for Ir, 3.92 A for Pt, 5.58 A for Ca, 
6.08 A for Sr, 4.05 A for Al and finally 4.95 A for P\xB. These numbers differ around 0.1 A from the numbers used 
in reference ^ where the theoretical GGA equilibrium lattice constants have been used. 

B. Stability of Anisotropy Ratios 

To test our convergence we present in table ^ the scalar-relativistic ASA low-index surface energies of Ag and the 
anisotropy ratios with respect to the different parameters used in the program. The absolute values of the energies 
change less than 0.01 eV and the anisotropy ratios change by less than 1%. The largest effect comes from the imax 
cut-off for the wavef unctions, but globally the first set of parameters is sufficient to give accurate values of both the 
surface energies and the anisotropy ratios. The second test presented in table || concerns the effect of the lattice 
parameter on the surface energies and on the anisotropy ratios. Our test has been also performed for the noble metals 
in scalar-relativistic ASA. Using Cu as a test case, the absolute values of the surface energies change by the same 
percentage for all the surface orientations when the lattice parameter is decreased. In the case of Cu the theoretical 
LDA lattice constant is around 2% smaller than the experimental one, so that this effect changes the surface energy 
by less than 2.5%. However due to error cancellation the anisotropy ratio changed by less than 1.1%. For Ag and in 
general for the other Ad metals, Rh and Pd, the LDA lattice constants are somewhat 1% smaller than the experimental 
lattice constants, while for the 5d metals the LDA values agree with the experimental lattice constants. Therefore the 
anisotropy ratios are independent of whether we use the LDA or the experimental lattice parameter and the following 
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discussion of the anisotropy ratios would not change if we would have used the LDA instead of the experimental 
lattice constants in our calculations. 



C. Relativistic Effects 

To study the importance of the relativistic effects we present in figure the non-relativistic (NR) and scalar 
relativistic (SR) density of states (DOS) for Au and Pt for the central layer atoms of the (111) slab which represent 
accurately the bulk DOS. Au has all d states filled and the c?-bands are deeper in energy compared to Pt that has only 
9 d-electrons and thus has a peak near the Fermi level. This means that the bonds between Pt atoms are stronger 
than for Au, so that for Pt more energy is needed to create a surface of a given orientation than for Au. Relativistic 
effects are small for the Pt DOS at the Fermi level but they widen the d-band so that the energy needed to break a 
bond is considerably increased. Although the SR-DOS of Au at the Fermi level is barely enhanced, relativistic effects 
broaden the band and shift it higher in energy thus enhancing the bonds between Au atoms and increasing the surface 
energies. Ag and Cu present a similar DOS with Au although the bandwidth of their d-bands is somehow smaller. In 
these cases, relativity only slightly changes the DOS, so the effect on the surface energies is smaller than for Au. 

In figure || we present the FP surface energies of the three noble metals obtained both in NR and SR calculations. 
The main relativistic effect is to increase the surface energies and this effect is largest for Au being the heaviest 
element. In the NR calculations Au has the lowest surface energies of all three noble metals but in the SR calculations 
Au sutface energies become comparable to the one of Cu and the inclusion of spin-orbit coupling further increases 
themcj. Also the anisotropy ratios are increased compared to non-rclativistic values, by about 2-4% in the case of Cu 
and Ag and up to 8-10% in the case of Au. 



III. SURFACE ENERGIES AND ANISOTROPY RATIOS 



We have studied the surface energies of the low-index surfaces, (111), (100) and (110) and of the three most close- 
packed vicinal surfaces: (311), (331) and (210). From a slab calculation one can calculate the surface excess free 
energy at zero temperature from the relation 

Eslab — NEbulk /-, s 

7= ^ (1) 

where Esiab is the total energy of the slab, N is the total number of layers of the metal, Ebuik is the per atom energy 
in the bulk crystal and the 2 enters because when we do a slab calculation we open two surfaces. To be consistent for 
all the cases we used as Ebuik t^e energy per atom of the central layer of the slab. 

The broken bond- rule statea£j that the surface energy ^(^hki) in eV/ (surface atom) needed to create a surface with 
a Miller index {hkl) reduces just to the product of 7(iii) and the ratio of the first-neighbor broken bonds iV(/jfe;) and 
iV(in) = 3: 

^ ^{hkl) 

lihki) = — ^ — 7(111)- (^) 

N^(hki) can be easily obtained for any fee surfaceEl: 

AT { 2h + k h,k, I odd l ; , fo\ 

N(hki) = [^h + 2k otherwise h>k>l. (3) 

As a consequence the anisotropy ratios between the surface energy for an arbitrary surface orientation over the 
surface energy for the (111) surface orientation is close to the ratio between the number of broken bonds between 
nearest neighbors for these surfaces. Here we will study how well the paramagnetic fee transition and sp metals satisfy 
the broken bond rule for the surface energies. 



A. Transition Metals 



In table III we have gathered our FP scalar relativistic results for all the transition metal surfaces. The first 



comment that someone can make on this table is pretty obvious. For all the metals the surface energy increases with 
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the roughness of the surface, i.e. as the number N^^/^f^^ of broken bonds increases, where (111) is the most close-packed 
surface with 7V(iii)=3. Also it is pretty obvious that surface energies for an isoelectronic row increase with the extent 
of the wavefunction; the surface energy of Pt is larger than the one of Pd, since the Pt 5d wavefunctions have a 
somewhat larger extent than the Pd 4c? ones. Also along a line surface energy is larger for the compound with the 
smaller number of d electrons as this one presents a stronger peak at the Fermi level (see figure^, e.g. Pt presents 
larger surface energies than Au. Also in the same table we present in parenthesis the anisotropy ratios. To open 
the (100) surface we break 4 nearest-neighbor bonds, the (110) 6 bonds, the (311) 7 bonds, the (331) 9 bonds and 
finally the (210) 10 bonds. So the ideal broken bond ratios with respect to the (111) surface, for which we break 3 
nearest-neighbors bonds, are 4/3, 2, 7/3, 3 and 10/3, for the (100), (HO), (311), (331) and (210) surface orientations 
respectively. The calculated surface energy anisotropy ratios deviate slightly from these ideal numbers. For Pd the 
ratios are smaller by the ideal ones by ~3-4% for all the surface orientations, while Ir ratios are larger from the ideals 
ones by ~4-7% for all the surface orientations except the (331) where the Ir calculated ratio is only by 1.3% larger 
than the ideal ratio. Pt and Rh show a mixed behavior but in general the ratios differ less than 3% from the ideal 
nearest-neighbors broken bonds ratios for any surface. So in general transition metal surfaces follow the broken bond 
rule but with slightly larger deviations than the noble metals due to the fact that their d band is not filled and they 
present peaks at the Fermi level, which can slightly change from one surface orientation to the other and consequently 
the energy needed to break a bond changes also slightly. 

In figure ^ we have plotted the surface energy anisotropy ratios with respect to the (111) surface for the low-index 
surface together with the results by Vitos et at from the reference In contrast to the noble metals the values of 
Vitos et al. only slightly deviate from our results for the transition metals. In reference |2^ we have explained the 
discrepancy for the noble metals as due to an insufficient number of k|| -points used in reference |2^ in the evaluation 
of the Brillouin zone integrals. The (111) surface of the noble metals presents a surface state centered at the f point 
which can only be accounted for by a sufficiently dense grid. Such states do not occur for the fee transition metals 
and thus the number of k|| -points used in reference 2q is sufficient to produce accurate ratios, that agree well with 
our FP results. Since the authors of the reference |25 did not perform calculations for the vicinal surfaces, we cannot 
judge the behavior of the FCD-LMTO-ASA method in these more difficult cases. 

As mentioned in the introduction there are several other ab-initio calculations for the surface energies of these 
materials but they concern just one or two surfaces and thus allow no conclusions for the anisotropy ratios. In table 
^ we have gathered the results from previous ab-initio calculations and experiments. We have expressed all the 
results in Jm~^ and not in eV/atom, as experiment have been performed in the liquid phase of the metals. For the 
noble metals our and Vitos' results agree nicely when expressed in Jm^^, and the calculated values are very close to 
the experimental ones. The latter values are not orientation specific but averaged values and they should be closer 
to the most close-packed surface: the (1-It1)- Other calculations agree with our and Vitos databases. The FP-LMTO 
results by Methfessel and coUaboratorsllj agree nicely with ours except the case of Ag where they predict a jellium 
like behavior for the (111) and (100) surfaces, i.e. same surface energy per unit surface area, which is not expected 
for a noble metal. Finally our [Calculated surface energies are in reasonable agreement with previous tight-bindiM 
calculations by Barreteau et alx3 on the low-index surfaces of Rh and Pd and by Mehl and Papaconstantopouloaij 
on the ensemble of noble and transition metals. 

We should mention that in our calculations we did not relax the positiopS-Df the layers but to a large extent 
this effect does not affect the surface energiefj. Feibelman and coUaboratorscHlEJ and Mansfield and collaborator^ 
showed by first-principle calculations that the effect of the relaxation on the calculated surface energy of a particular 
facet should be around 2-5% depending on the roughness of the facet. Surface relaxations for vicinal surfaces have 
been studied mainly usiag. semi-empirical methods due to the complexity arisen by the simultaneous relaxation of 
a large number of layersE^. Rodriguez and collaborators using jSuch a semi-empirical method showed that surface 
relaxations affect typically the anisotropic ratios by less than 2 %Ej, and so the neglect of relaxation should have little 
effect on our results. 



1. ASA versus FP Calculations 



The results presented in the previous section have been obtained taking into account the non-spherical part of the 
potential, i.e. known as the full-potential scheme. The use of the FP instead of the ASA accounts in a more accurate 
way for the charge distribution near the surface where due to the lower symmetry the charge exhibits larger variations 
than in the bulk. The use of FP lowers the energies compared to ASA by about 15% for all the surfaces under study. 
A similar behavior is also found for the vicinal surfaces. At first sight it seems that the ASA is efficiently accurate to 
describe the surface energies of the materials under study. However, the anisotropy ratios are more sensitive that the 
surface energies themselves. 



4 



In figure ^ we have represented the anisotropy ratios for the noble and transition metals for the three low index 
and the three vicinal surfaces which we have studied. The straight lines represent the ideal broken bonds ratios. We 
see that already for the Pd(311) surface the ASA produces a ratio that deviates strongly from the broken bond rule 
while FP is near it. The differences become even more dramatic for the more open (331) and (210) surfaces. ASA 
produces ratios for the Cu and Ir surfaces that deviate strongly from the broken bond rule while using the FP we find 
again the ideal ratios. Especially for Pd, the ASA predicts that the surface energy for the (331) surface is larger than 
the one for the (210) surface which is more open. The use of FP restores the broken-bond rule behavior. So although 
the ASA is sufficiently accurate to produce reasonable surface energies, the use of FP is decisive for the calculation 
of the anisotropy ratios of the vicinal surfaces. 



B. sp Metals 



To complete our study we also investigated the paramagnetic sp metals - Ca, Sr, Al and Pb - that crystallize also 
in the fee structure. In table ^ we have gathered the FP surface energies and in parenthesis the anisotropy ratios for 
the six more close-packed surfaces. In general the surface energies for these materials are smaller than for the d metals 
due to the fact that the bonds are made of s and p electrons that are more mobile than the localized d electrons and 
so one needs less energy to break these bonds. This becomes even more clear when we look at the surface energy 
expressed in surface units (see table [v| ). The Ca(lOO) and Sr(lOO) surfaces show a rather small anisotropy ratio 
compared to the transition and noble metals and also compared to Al. Ca and Sr. are in the periodic table just near 
the simple metals which are known to be well described by a simple jellium modcQ. So we expect that to some extent 
we would find a jellium like behavior also for the Ca and Sr surfaces at least for the low-index ones. This is really 
what happens. In figure ^ we represent the anisotropy ratios but now taking into account the energy per surface area 
(Jm~^) and not per atom (eV/atom). For a jellium model the surface energy per surface area would be constant for 
any surface and the anisotropy ratio would be always 1. With the solid line we represent the anisotropy ratios if the 
broken-bond rule is applicable, and the ratios for Ca, Ag and Ir. Ag is closer to the broken bond rule than Ir where 
the ratios are always overestimated. For the Ca(lOO) surface we see that the value is closer to the jellium model but 
for the more open surfaces the anisotropy ratios are closer to the broken bond model than the jellium. 

Contrary to Ca and Sr, Al shows the same behavior with Pd and the calculated ratios are slightly smaller than the 
ideal ones (see table 0). The most interesting case is Pb. For the (100) and (110) surfaces the ratios are near the 
Sr ones while for the next two vicinal surface orientations, (311) and (331), the calculated anisotropy ratios are more 
than 6% smaller than the ideal values. But for the (210) surface the anisotropy ratio changes the behavior and now is 
larger by 6% than the ideal value of 10/3. This is an indication that the behavior of Pb is more complicated than all 
the other fee metals we have studied in this contribution. Here we have to mention that the surface energies in table 
have been calculated taking the 5d as valence electrons as they are located just below the sp bands. We recalculated 
the surface energies of Pb considering the 5d as core states, by increasing the number of layers, the number of energy 
points used to do the integrations in the complex energy plane and finally the number of k||-points. Although the 
surface energies changed slightly due to the core treatment of the 5d states, the anisotropy ratios were extremely 
stable. 

In table VI we have gathered our calculated surface energies in Jm^^ together with other calculations and exper- 
iments. Our results agree nicely with the Vitos et al. database in reference |2^ with the exception of Pb where our 
surface enemy per surface area is double as high, but our results agree nicely with previous calculations by Mansfield 
and NeedsLJ using pseudopotentials and the existing experimental data. Both Schochlin et a/.cEl and Stumpf and 
Schefflercj have studied all the three low-index surfaces of Al. Comparing Schochlin et al. calculations with Stumpf's 
and Scheffler's results we see that the former calculations predict comparable surface energies for the (100) and (110) 
surfaces while the latter ones predict comparable surface energies for the (111) and (100) surfaces. Both calculations 
are in contradiction to our calculations that predict a considerable increase of the surface energy as the surface be- 
comes more open, while the calculations in reference |2^ state that the surface energy per surface area is smaller for 
the (110) surface compared to the (100) surface. This spread of the results for Al does not allow us to draw safe 
conclusions for the variation of the surface energy with the surface orientation. 



IV. CONCLUSIONS 



We have shown using the full-potential (FP) screened KKR code that the broken bond rule, i.e. the surface energy 
scales linearly with the number of nearest-neighbors broken bonds, already shown for the noble metals is also valid 
for the transition metals. For the sp metals the (111) and (100) surfaces show a jellium like behavior but for the more 
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open surfaces the surface energies follow again the broken-bond rule with Pb presenting the largest deviations. The 
use of full-potential instead of the full-charge atomic-sphere approximation (ASA) decreased the surface energies and 
we showed that it is necessary to accurately calculated the anisotropy ratios for the vicinal surfaces. 
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TABLE I. Scalax-relativistic surface energies and anisotropy ratios in parenthesis for Ag within ASA using different values 
for the the imax cut-oil (LM), the tight-binding cluster (TB), the number of energy points (EN) to perform integrations in the 
complex energy plane and the number of k||-points in the full two-dimensional Brillouin zone. The first set of parameters is 
the one used in the present calculations. 

LM TB EN KP Ag(lll) Ag(lOO) Ag(llO) 

~3 19 27 55 x 55 0.641 (1.34) 0.860 (1.98) 1.271 

3 55 27 55x55 0.637 (1.34) 0.854 (1.98) 1.262 

3 19 41 55x55 0.643 (1.33) 0.855 (1.97) 1.266 

3 19 27 75x75 0.648 (1.33) 0.860 (1.96) 1.273 

3 55 41 75x75 0.642 (1.32) 0.849 (1.96) 1.259 

4 19 27 55x55 0.649 (1.35) 0.879 (2.00) 1.301 



TABLE IL Effect of lattice parameter changes on the scalar-relativistic ASA surface energies for the three low-index surfaces 
of Cu in the upper panel and for the (111) surface of Cu, Ag and Au in the bottom panel, 
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TABLE in. Pull-potential scalar-relativistic surface energies for the six more closed-packed surfaces for the four transition 
metals. In parenthesis the anisotropy ratios with respect to the (111) surface. 



7(eV) Pd Pt Rh Ir 

(111) 0.822 0.957 1.034 1.200 

(100) (1.28) 1.049 (1.33) 1.272 (1.36) 1.404 (1.42) 1.707 

(110) (1.94) 1.596 (2.06) 1.973 (1.98) 2.047 (2.07)2.488 

(113) (2.28) 1.873 (2.40) 2.295 (2.35) 2.428 (2.43) 2.913 

(331) (2.93) 2.404 (2.98) 2.853 (2.99) 3.094 (3.04) 3.652 

(210) (3.22) 2.644 (3.30) 3.158 (3.35) 3.464 (3.48) 4.172 
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TABLE IV. Surface energies is J m ^ for the three low-index surfaces of the paramagnetic fee transition metals using the 



FKKR. In second column the results using the FCD-LMTO-ASA and on the third and fourth columns existing ab-mitio 
calculations and experiments. 
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TABLE V. Full-potential scalar-relativistic surface energies for the six more closed-packed surfaces of the four sp metal. In 
parenthesis the anisotropy ratios with respect to the (111) surface. 



7(eV) Ca Sr Al Pb 

{111) 0.417 0.373 0.489 0.398 

(100) (1.21) 0.503 (1.22) 0.454 (1.28) 0.625 (1.24) 0.492 

(110) (1.92) 0.797 (1.94) 0.722 (1.93) 0.943 (1.95) 0.778 

(311) (2.27) 0.946 (2.27) 0.847 (2.24) 1.094 (2.18) 0.866 

(331) (2.96) 1.234 (2.93) 1.094 (2.94) 1-436 (2.81) 1.118 

(210) (3.20) 1.333 (3.18) 1.187 (3.20) 1.565 (3.53) 1.405 
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TABLE VI. Surface energies is Jm~^ for the three low-index surfaces of the paramagnetic fee sp metals using the FKKR. 
In second column the results using the FCD-LMTO-ASA and on the third and fourth columns existing ab-imtio calculations 
and experiments. 
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FIG. 1. Non-relativistic (NR) density of states (DOS) of Pt (upper panel) and Au (bottom panel) compared with the 
scalar-relativistic (SR) DOS. All DOSs are calculated for the central layer atom of the (111) slab. 
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FIG. 2. Effect of relativity on the full-potential surface energies of the noble metals. The effect of relativity is large for Au 
that possesses 5d electrons. Relativity also enhances the anisotropy ratios. 
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FIG. 3. Anisotropy ratios within both FKKR and LMTO from reference 25. Contrary to the noble metals, whose (111) 
surface posses a surface state centered at the F point, the transition metals results are similar for both methods. 
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FIG. 4. Anisotropy ratios for the low-index and the most close-packed vicinal surfaces within ASA and FP. For the vicinal 
surfaces the use of FP is necessary, due to the failure of ASA to describe accurately the more complex metal- vacuum interface. 
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FIG. 5. Anisotropy ratios but calculating the surface energies in Jm~^ for Ag, Ir and Ca. The solid line corresponds to the 
brokon-bon rule and the tilted one to the jellium. The Ca(lOO) surface is close to jellium but then it recovers the Ag and Ir 
behavior. 
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